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CHARTS  FOR  THE  PARAMETERS  OF  MIGRATING  EXPLOSION  HUBBIES 

By 


H.  G.  Snay  and  Ruth  V.  Tipcon 


ABSTRACT:  Graphs  are  presented  for  the  relative  bubble  energies,  maximum 
radii,  periods,  and  migrations  of  underwater  explosion  bubbles  for  four 
cycles  of  the  bubble  oscillation.  The  migration  refers  to  the  bubble 
maximum  in  each  cycle.  The  position  of  the  bubble  minimum  can  be  readily 
found  by  a  simple  confutation.  The  graphs  permit  a  convenient  reading  of 
these  parameters  for  a  wide  range  of  conditions. 
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CHARTS  FOR  TOE  PARAMETERS  OF  MIGRATING  EXPLOSION  BUBBLES 

This  report  is  part  of  a  comprehensive  study  on  the  behavior  of  explosion 
gas  bubbles  which  has  been  carried  out  tinder  Task  No.  30I-66VU3003 /OlO^O 
and  HUME-3-E-000/212-l/WF  008-10-004  PA  002.  A  knowledge  of  the  behavior 
of  migrating  gas  bubbles  is  of  considerable  importance  in  studies  of  ship 
damage.  NAVORD  Report  3906,  "Key  Problems  in  Explosive  Research  and 
Development",  Chapter  III,  "Anti-Submarine  Warfare"  discusses  this  in  more 
detail.  The  present  report  represents  a  numerical  evaluation  of  the  theory 
developed  in  NAVORD  Report  4l85.  The  charts  prepared  on  the  basis  of  this 
theory  provide  a  simple  method  which  this  Laboratory  believes  will  save 
considerable  time  and  effort  in  calculating  the  important  characteristics 
of  explosion  bubbles  in  the  second  and  subsequent  cycles  of  the  pulsation. 

The  major  portion  of  this  work  was  carried  out  during  a  course  of  duty  of 
the  second  author  as  a  Junior  professional  trainee. 

R.  E.  ODENING 
Captain,  USN 
Commander 

C.  J.  ARONSON 
By  direction 
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I.  INTRODUCTION 

In  NAVORD  Report  4l85*  relationships  are  derived  for  the  migration  and 
energy  ratios  of  underwater  explosion  bubbles.  The  data  given  in  that  re¬ 
port  are  sufficient  to  provide  the  information  needed  for  calculating  the 
migration  period,  maximum  radius,  and  energy  through  the  4th  cycle.  In  the 
present  report  this  information  has  been  used  to  make  graphs  for  the  ratios 
of  these  parameters.  The  use  of  non-dimensional  variables  permits  the 
convenient  reading  of  the  results  for  a  wide  range  of  explosive  conditions. 

II.  THE  EQUATIONS  OP  MIGRATING  BUBBLES 

The  equation  for  the  migration  for  pulsating  bubbles  between  the  nth 
and  (n  +  l)th  bubble  maximum  derived  in  NAVORD  Report  4l85  (Equation  9)  can 
be  written  in  the  following  form 


Z  =■  Total  hydrostatic  pressure  measured  in  feet  of  water. 

1  atmosphere  »  33  ft  sea  water  =  34  ft  fresh  water 

^  *  Maximum  bubble  radius  in  ft. 

r  a  Ratio  of  bubble  energy  to  the  total  energy  of  explosion, 

n  ■»  Subscript  designating  the  cycle  of  the  pulsation. 

D  =  Depth  in  ft  below  the  water  surface.  D  refers  to  the  center 
of  the  bubble  maximum.  In  the  first  cycle.  D-,  1b  assumed  to 
be  equal  to  the  depth  of  explosion,  i.e.  the  migration  up  to 
the  first  bubble  maximum  is  neglected. 


*  Full  reference  is  given  at  the  end  of  this  report. 
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It  is  stated  in  NAVQRD  4l85  that  C  can  he  considered  to  he  a  universal 
constant,  C  =  3*5- 


From(l)  it  appears  that  LZ^/Z-^  is  a  function  of  Aj^/Z^,  ZJ1/Z1,  Dn/Zi 
and  r  /r. .  The  energy  ratio  rn/r-^  is  evaluated  in  NAVORD  4l85  and  is  shown 
to  he  a  function  of  AZn/Zn.  The  two  variables  2^/Z^  and  Dn/Z1  are  functions 
of  the  reduced  migration  AZq/Z^,  as  can  he  seen  in  the  following  relations: 


and 


gn  _  Zn-1  . 

Z1  Z1  Z1 

n  »  2,  3,  ... 

z,  Zn  Z, 


(2) 

(3) 


Therefore,  the  reduced  migration  and  the  depths  of  the  later  huhhle  maxima 
are  functions  of  two  variables  only,  namely  A^/Z-j^  and  D. /Z..  This  fact 
permits  a  convenient  graphical  representation  not  only  or  the  migration,  hut 
also  of  the  other  huhble  parameters. 


The  equation  for  the  first  maximum  huhhle  radius  is,  reference  (b), 


/  W  n 

AjO.  -  J  V  Zi  ) 


where 

W  *  Charge  weight  in  lbs, 

J  =  Radius  constant  dependent  on  the  properties  of  the  explosive, 
e.g.  J  =  12.6  for  TOT,  reference  (c). 

The  maximum  radius  at  any  one  of  the  subsequent  cycles  is  given  by 

Sfa  (  r_n 
^  "  Vrl 

The  period  of  the  first  bubble  pulsation  for  an  explosion  in  infinitely 
deep  water  is,  reference  (b), 


Z  A/3 

_1  ^  ,  (NAVORD  4185,  Equation  4). 

Z  J  V? > 
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wl/3 

7375 


i  - 


O.l 


^Ml  "I 


(6) 


where  K  is  the  period  constant,  e.g.  K  =  4.36  for  TNT,  reference  (c).  The 
term  0.1  Aj^,/'Dn  accounts  for  the  effect  of  the  water  surface  on  the  pulsa¬ 
tion.  It  is  negligibly  small  for  a  deep  explosion,  hut  compensates  for  the 
shorter  periods  of  a  bubble  pulsating  near  the  water  surface.  (Compare 
Section  III  of  NAVORD  4l85  about  the  value  0.1.) 


The  periods  of  the  subsequent  cycles  can  be  obtained  from 


T  _  l/3  z  ^5/6  1  -  0.1  A,,  /D 

rn  V  f  1  Y  _ Mn  a 

V  rj  J  K  2n  7  1-0.1  A^/Di 


(7) 


as  discussed  in  NAVORD  4185,  Equation  (2). 


Since  (5)  and  (7)  contain  only  parameters  which  are  functions  of 
Afto/Zi  and  D./Z^,  the  period  and  radius  ratios  are  functions  of  these  two 
variables  only  and  each  can  be  graphically  presented  by  means  of  a  family 
of  curves. 


For  shallow  explosions,  the  presence  of  the  water  surface  perturbs  the 
bubble  oscillation  in  such  a  way  that  the  emission  of  the  bubble  pulse 
ceases.  This  phenomenon  is  commonly  called  venting  of  the  bubble  or  blow- 
in.  It  can  be  quantitatively  taken  Into  account  using  the  evidence  that 
the  nth  bubble  pulse  will  not  occur  when 


Aj*m  s  9b* 


(8) 


III.  NUMERICAL  CALCULATIONS 

The  calculations  of  this  paper  are  based  on  the  relationships  listed 
above  and  the  values  of  r  / r,  found  in  NAVORD  Report  4l85.  The  ratios  for 
the  depths  of  the  bubble  maximum,  the  periods,  the  maximum  radii,  and  the 
bubble  energies  were  calculated  for  the  first  four  cycles  of  the  oscilla¬ 
tion.  In  view  of  the  uncertainty  of  the  energy  ratio  for  the  5th  cycle, 
given  in  NAVORD  4185,  this  case  was  not  included.  Graphs  were  made  using 
Aj^/Z,  as  abscissas  and  the  above  bubble  parameters  as  ordinates.  The 
result  is  presented  as  a  family  of  curves  each  labelled  by  a  value  D 
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Limiting  curves  were  constructed  for  the  condition  where  hlow-in  occurs, 
Figures  1  to  3.  In  the  other  charts  the  points  of  hlow-in  are  marked.  The 
continuation  of  the  curves  beyond  these  points  is  shown  in  dashed  lines. 
They  are  useful  for  certain  calculations,  although  for  such  conditions  no 
bubble  pulse  occurs. 


IV.  USE  OF  THE  GRAPHS 


Given  the  weight  of  the  explosive  charge  W  and  the  depth  at  which  it 
is  fired,  one  uses  the  formulae 


21 


ftj_  +  33 


(9) 


ha 

Zn 


-  J 


wl/3 

ZV3 


(10) 


to  find  the  input  of  the  charts.  The  ratios  read  from  the  graphs  are 
multiplied  with  the  value  of  the  corresponding  variable  for  the  first  cycle 
of  the  oscillation.  This  immediately  yields  the  desired  result. 

As  an  example  suppose  one  desires  to  know  the  depth  of  the  third  bubble 
maximum  of  a  290-lb  TNT  charge  fired  at  a  depth  of  150  ft  in  sea  water. 

Here,  =  150  ft  +  33  ft  =  183  ft.  Using  (10)  one  finds  A^j/Z.  m  0.0803 

and  from  (9)  Dj/Z^  =  0.82.  The  first  period  is  according  to  (6)  T,  =  O.372. 
Using  these  values  one  finds  from  Figure  2:  D^/D.  =  0.030  and  Do  =  125  ft. 

Thus,  the  upward  migration  of  the  bubble  from  the  initial  charge  location 

up  to  the  moment  of  the  third  bubble  maximum  has  been  25  ft.  In  an  analogous 

way  one  finds  from  Figure  5:  To/T^  =  . 75 >  and  To  =  0.2T9  seconds.  The  period 
T3  refers  to  the  time  between  the  second  and  third  bubble  minimum.  The  time 
or  the  third  maximum  would  be  T^  +  Tg  +  Tj/2  =  0.83*1-  seconds.  The  position 
of  the  third  bubble  minimum  can  be  found  By  calculating  that  of  the  fourth 
bubble  maximum  and  placing  the  bubble  minimum  half  way  between  these  two 
maxima. 


V.  COMMENTS  ON  THE  DEPENDENCY  OF  THE  BUBBLE  PARAMETERS 
ON  THE  EFFECT  OF  THE  WATER  SURFACE 

Inspection  of  the  graphs  presented  in  this  report  gives  an  interesting 
insight  on  the  influence  which  the  two  input  variables  have  on  the  bubble 
behavior. 
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The  magnitude  Aj^/Z^  Is  a  measure  of  the  buoyancy  of  the  bubble.  The 
maximum  radius  used  in  this  ratio  represents  the  pressure  difference  between 
bubble  top  and  center.  It  is  this  type  of  pressure  difference  which,  when 
integrated  over  the  bubble  surface,  produces  buoyancy.  stands  for  the 
total  pressure  at  the  bubble  center.  Hence,  Aj^/Z^  represents  one-half  of 
the  pressure  difference  between  bubble  top  and  bottom  divided  by  the  pressure 
at  the  center.  This  magnitude  accounts  for  the  driving  force  of  the  gravity 
migration. 

The  other  variable,  D^/Zp  is  a  measure  of  the  depth  at  which  the  explo¬ 
sion  occurs.  The  deeper  the  explosion, the  closer  is  this  value  to  unity. 

Comparing  the  above  equations,  it  is  seen  that  D^/Z^  occurs  in  the 
surface  correction  terms  only.  Hence,  Aj^/Z^  is  a  measure  of  the  strength 
of  migration,  and  D^/Z^  is  a  measure  of  the  surface  effect. 

The  graphs  illustrate  the  importance  of  these  two  effects  on  the  bubble 
parameters.  Figures  1  to  3  represent  the  migration.  The  large  spread  of  the 
family  of  curves  in  this  figure  indicates  that  the  surface  effects  have  a 
strong  bearing  on  migration.  The  influence  is  less  pronounced  for  the 
bubble  periods,  Figures  4  to  6,  but  still  clearly  noticeable.  This  evidence 
was  to  be  expected,  since  equations  (l)  and  (6)  include  surface  correction 
terms.  But,  the  surface  effect  also  enters  into  the  computation  of  the 
maximum  bubble  radius  and  the  bubble  energy,  although  these  terms  are  not 
directly  affected  by  the  water  surface.  The  curves  for  the  maximum  radii, 
Figure  7,  show  this  influence  only  for  close  proximities  of  the  bubble  to 
the  water  surface:  The  curve  for  each  reduced  depth  D. /z  emerges  from  the 
average  curve  near  the  point  of  blow-in  or  for  very  shallow  conditions.  The 
same  effect  is  noticeable  for  the  bubble  energy  ratio,  Figure  8. 

This  discussion  illustrates  the  importance  of  the  surface  effect  and 
shows  that  it  must  not  be  neglected  in  a  study  of  bubble  migration.  The 
charts  presented  in  this  paper  quantitatively  account  for  this  effect  and, 
at  the  same  time,  allow  for  a  convenient  appraisal  of  its  significance.  For 
instance  maximum  radius,  bubble  energy,  and  periods  depend  mainly  on  Aj^/Zp 
On  this  basis,  simple  approximate  formulae  could  be  set  up  which  give  these 
magnitudes  as  functions  of  Aj^/Z^  only.  The  charts  illustrate  the  limited 
accuracy  of  such  approximations.  Since  reading  of  a  graph  is  simpler  than 
to  evaluate  equations,  we  have  not  included  such  formulae  in  this  report. 
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FIG.  I  RELATIVE  DEPTH  AT  SECOND  BUBBLE  MAXIMUM 
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FIG.  2  RELATIVE  DEPTH  AT  THIRD  BUBBLE  MAXIMUM 
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FIG.  4  RELATIVE  SECOND  BUBBLE  PERIOD 
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FIG.  5  RELATIVE  THIRD  BUBBLE  PERIOD 
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FIG.  6  RELATIVE  FOURTH  BUBBLE  PERIOD 
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RATIO  OF  MAXIMUM  BUBBLE  RADII 
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FIG.  7  RELATIVE  MAXIMUM  BUBBLE  RADII 

D,/Z,  LABELS  ARE  ONLY  SHOWN  FOR  THE 
UPPER  FAMILY  OF  CURVES.  LABELS  FOR  THE 
OTHER  CURVES  ARE  ANALOGOUS 
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RATIO  OF  BUBBLE  ENERGIES 
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FIG,  8  RELATIVE  BUBBLE  ENERGIES 
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